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Abstract 

Helically  syinitietric  hydromaynetic  equilibria  with  small 
helical  corrugation  are  considered  in  a  finite  cylindrical 
domain.   These  solutions  may  be  regarded  as  bifurcated  from  a 
circularly  symmetric  equilibrium  with  the  plasma  current  at  the 
plasma  surface.   Stable  bifurcated  equilibria  must  have  a 
bifurcation  point  in  the  neighborhood  of  the  margin  of  stability 
for  the  bifurcating  equilibria.   Using  the  energy  principle  we 
show  that  there  are   stable  helical  equilibria  with  long  helical 
wavelength,  the  corresponding  circular  equilibria  are  unstable. 
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I .   Introduction 

Bifurcation  of  solutions  representing  equilibria  occurs  in 
many  physical  problems.    A  well-known  example  is  the  buckling 
of  an  elastic  beam.   As  the  characterizing  parameters  vary, 
one  family  of  solutions  (the  bifurcating  equilibria)  changes 
from  stable  to  unstable,  and  the  other  family  of  solutions  (the 
bifurcated  equilibria)  branches  out  at  the  bifurcation  point. 
The  bifurcated  equilibria  are  usually  of  lesser  symmetry  than 
the  bifurcating  equilibria.   Since  solutions  representing  stable 
hydromagnetic  equilibria  are  important  in  magnetic  confinement 
of  a  plasma,  it  is  of  interest  to  study  the  stability  of 
bifurcated  hydromagnetic  equilibria. 

The  simplest  hydromagnetic  model  for  a  closed  system  of  plasma 
containment  is  a  straight  sheet  pinch.   Namely,  inside  a 
circular  cylinder  of  conducting  wall  the  plasma  forms  a  column 
Vv'hose  cross-section  is  a  circle  concentric  with  the  outer  wall  . 
In  the  plasma  region,  the  magnetic  field  has  only  an  axial 
component,  and  in  the  vacuum  region  between  the  plasma  and  the 
conducting  wall  the  magnetic  field  has  both  an  azimuthal  as  well 
as  an  axial  component.   Besides  the  size  (as;  specified  by  its 
length  and  radius)  of  the  wall,  under  the  assumption  of  a  sharp 
boundary  for  the  plasma-vacuum  interface   the  system  is  character- 
ized by  four  physical  parameters:  the  axial  magnetic  flux  in  the 
plasma,  the  axial  flux  in  the  vacuum,  the  total  azimuthal  flux, 
and  the  volume  of  the  plasma.   Given  the  four  parameters,  in 
addition  to  the  circular  equilibrium  helical  equilibria  are 
other  possible  configurations.   The  helical  structure  is  specified 
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by  an  axial  wavenumber ,  an  azimuthal  wavenumber,  and  the 
amplitude  of  helical  corrugation.   The  admissible  axial 
wavenumbers  are  determined  by  the  length  of  the  conducting  wall. 
The  amplitude  of  corrugation  is  determined  by  the  physical 
parameters . 

The  bifurcation  problem  for  a  hydromagnetic  system  v/as  first 

2 
discussed  by  Friedrichs.    He  considered  the  same  helical  equi- 
libria as  is   discussed  in  this  paper   in  the  limit  of  large 
wavenumbers  (both  axial  and  azimuthal).   Clearly,  the  helical 
equilibria  with  small  amplitude  can  be  stable  only  when  the 
bifurcation  point  is  in  the  neighborhood  of  the  margin  of 
stability  for  the  corresponding  circular  equilibria.   This 
requires  the  azimuthal  wavenumber   to  be  equal  to  unity.   Also 
it  is  clear  that  any  instability  for  such  bifurcated  equilibria 
must  arise  from  disturbances  with  an  azimuthal  wavenumber  equal 
to  unity  and  an  axial  wavenumber  differing  at  most  slightly 
from  that  of  the  helical  structure. 

In  Section  II  we  consider  helically  symmetric  equilibria, 
expanding  the  solution  in  a  small  parameter,  the  amplitude  of 
helical  corrugation  of  the  plasma-vacuum  interface.   In 
Section  III  we  discuss  bifurcation  points  in  a  parameter  space. 

Then,  we  study  the  stability  of  the  bifurcated  helical  equilibria, 

3 
using  the  energy  principle  for  hydromagnetic  stability.   In 

minimizing  the  second  variation  of  the  potential  energy,  the 

displacement  is  compressible  if  the  disturbance  is  in  the  same 

mode  as  that  of  the  helical  structure,  and  is  incomprt-'ssible 

otherwise.   We  treat  the  two  cases  separately  in  Sections  IV  and  V, 


-2- 


Finally  in  Section  VI  we  show  that  there  are  stable  helical 
equilibria  with  long  helical  wavelength. 

Usually  when  the  characterizing  parameters  are  such  that  the 
circular  equilibrium  is  unstable, the  condition  is  interpreted  as 
unstable.   However,  if  there  exists  a  stable  helical  equilibrium 
for  the  same  parameters ,  then  the  correct  interpretation  should 
be  that  the  condition  is  stable  with  the  plasma  appearing  in  a 
helical  configuration. 

Recently,  hydromagnetic  equilibria  with  helical  symmetry 

4 
were  calculated  numerically  by  Friednan.    And  stability  for 

certain  equilibria  with  long  helical  wavelength  was  analyzed 

5 

by  Weitzner. 
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II.   Equilibria  with  Helical  Symmetry 

We  consider  a  helically  symmetric,  hydromagnetic  equilibrium, 
in  which  a  plasma  is  confined  by  magnetic  fields  inside  a 
circular  cylinder  of  conducting  wall.   The  plasma  current  is 
assumed  to  be  concentrated  at  the  surface  and  the  sharp  plasrfia- 
vacuum  interface  is  assumed  to  have  a  small  amplitude  of  helical 
corrugation. 

In  cylindrical  coordinates  (r,9,z)   the  helical  symme1-ry 
requires  that  all  quantities  depend  only  on  the  radial  variable 
r  and  the  helical  variable 

4)  =  m6  +  kz  , 

where  m  and  k  are  the  azimuthal  and  axial  wavenumbers  character- 
izing the  helical  structure  (m  can  be  any  integer  but  k  must  be 

one  of  the  discrete  values  determined  by  the  length  of  the  wall) . 

->• 

The  magnetic  field   R  =  (b  ,B„  ,B  )   can  be.  described  in  terms 

r  '  0  '  z 

-> 

of  a  flux  function  *i'(r ,({)).   The  requirement  V 'B  =  0  is 
satisfied  by  the  specification  B   =  -r~  9f/8(})  and 
riBg  +  krB^  =  Bf/ar.   Since   B«Vf  =  0,   '1'  is  constant  on  each 
flux  surface.   The  assumption  Vxb  =  0  requires  that  "F  satisfy 

,-|x  13  r  3T         1      3^4'         _  2mkA 

r    9r  ~T~iri:  3F  "^   "?  77T"     ~      .    2^,2    2,2    ' 
m  +k   r  r      3(|)  (m  +k   r    ) 

where    A    is    a   constant   equal    to   mB      -   krB„ .      Accordingly, 

(2)  ^    _    /       1    3*         "^    3F   -   ^^^         ^^    3F  +   ™M 

\~  r   W   '    2^,  2   2      '    2^^2    2      )    • 

^  m  +k   r  m  +k   r        ' 


Thus,  on  a  cylindrical  surface  r  =  const.   the  magnetic  field 

2   2  2—1  /2 
has  a  component  equal  to  Cm  +k  r  )   '^  A  along  the  helix 
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2   2  2  -1/2 

m9  +  kz  =  const,  and  a  component  equal  to  (m  +k  r  )  '     3f/8r 

perpendicular  to  the  helix.   The  boundary  conditions  to  be 
satisfied  are  that  the  plasma-vacuum  interface  and  the 
conducting  wall  are  two  flux  surfaces  and  the  jump  in  the 
magnetic  pressure  across  the  interface  is  balanced  by  the 
constant  plasma  pressure. 

The  helical  configuration  is  characterized  by  four  physical 
quantities: 
the  axial  magnetic  flux  in  the  plasma  region 


2u 


(3) 


f 


9=0   r=0 


B  r  dr  de  , 
z         ' 


the  axial  magnetic  flux  in  the  vacuum  region 


(4) 


B   E 


27T 

r 


e 

r 


B  r  dr  d( 
z 


)  =  0   r=p 


the  total  azimuthal  magnetic  flux 


(5) 


e 

r 


z=0   r=0 


B„  dr  dz 


and   the   volume   of    the   plasma 


2tt 


(6) 


z=o     e 


0 


J  p      d9    dz    , 


where  e  and  L  are  the  radius  and  the  length  of  the  cylindrical 
wall,  and  p  is  the  plasma-vacuum  interface:  r  -  p((J))  =  0. 

We  shall  use  the  following  notations  defined  in  terms  of  the 
standard  modified  Bessel  functions  of  integral  order   I  (x) , 
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K  (x)  and  their  derivatives   I  (x)  ,  K  (x)  : 

V  V        V 


1    ^v^^^^ 
I(r;v,K)  =  -^     — i , 


I^«d) 
T    K^(Kr)l'(Ke)  -  I,  (Kr)K'(Ke) 


K^(Kd)I^(Ke)  -  I^(Kd)K^(Ke) 


W{v,k)   e   (Kad)^  T(d;v,K)  -  c^  -  (vc+Kbd)^  K(d;v,<)  , 


I^^(r)  =     I(r;vm,vk)  ,  ^  o(b^^^     =  K(r;vm,vk)  ,   v  =  1,2,...  . 

In  expressions  which  are  clear  from  the  context  involving  no 
coordinate  variables,    I,  ,  ^  2<i)    '    ^  <h    '    ^  2(b      ^^^^  ^^   understood 
as  evaluated  at  r  =  d.   Also  we  shall  use  equal  signs  to  mean 
"equal  to  the  real  part  of"  whenever  the  right  side  is  complex- 
valued,  and  use  asterisks  to  mean  complex  conjugates. 

The  plasma-vacuum  interface  is  a  free  surface  with  helical 
structure.   When  the  helical  corrugation  is  small  p{<p)    can  be 
written 

p  =  d[l  +  e  exp  (i(}))  +  1  e^A  exp  (i2({))  +  0(e"^)  I  . 

The  small  parameter  e   is  a  measure  of  the  helical  corrugation, 

and  the  constant   X  is  to  be  determined  by  the  boundary  conditions _ 

In  the  plasma  region  (0  <  r  <  p)  we  take  A  =  ma  and  write  a  solution 
of  Eq.  (1)  as 


r^  d 

^  =  aQ  +  ka  2-  +  e  a^r  g^  I^  exp  {i<^) 

^  T  ^  ^2^  dr   ^2(J)  ^^P  ^^^*^  ^  °^^  ) 
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The  corresponding  magnetic  field  is 


.  d 


m 


(7) 


B  =  (0,0, a)  +  ea^(-i  ^     ^^    ^    ^     ^  <^    '    ^      ^  ^^    ^^P  ^^^^ 


+  J   ^S(-i  aF  ^2c^  '  ^  ^2   '  2k  l2^)  exp  (i2c^)  +  OCe^) 


In  the  vacuum  region  (p  <  r  <  e)  we  take  A  =  mb  -  kcd  and 
write  another  solution  of  Eq.  (1)  as 


2   2 

^  =  b^  +  mcd  log  |  +  kb  ^  ^^   +  eb^r  ^  K    exp  (i({)) 


0 


+  ^  £  b2r  ^  K^,  exp  (i2()))  +  0(e  )  . 


The  corresponding  magnetic  field  is 


B  =  (0,c  I  ,b)  +  eb^(-i  ^K       ,  ^K   k  K    )    exp  (i(j)) 


(8) 


-^  7   ^S(-i  aF  S*  '  ^  Scj)  '  2k  ^2^)  exp  (i2^)  +  O(e^) 


The  various  coefficients  will  be  determined  in  terms  of  the  four 
physical  parameters  A,  B,  C,  and  D  by  the  boundary  conditions. 


The  constants   a„  and  b-   are  related  to  the  magnetic  fluxes 


by 


^0  -  -  ^7 


•^0  =  -YW 


'■n  p 

T 

0     0 

2tt     e 

r 


(mB„  +  krB  )  dr  , 
9       z      ' 


d(j)  1  (mB   +  krB^)  dr 


From  the  requirement  that  ^  is  constant  at  the  interface  we  find 
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a-j^  =  -  kad^  +  O(e^)  , 


b,  =  -  (mc+kbd)d  +  0(e  )  , 


a^   =   -   kad' 


1  -  2(m^+k^d^)  I   +  2A 


+  0(e)  , 


h^   =  2mcd  -  (mc+kbd)d 


1  -  2(m^+k^d^)K   +  2x1  +  0(e) 


1  -»■  ^ 
2-B-B 


Next,  from  the  requirement  that 

equal  to  the  plasma  pressure  p,  we  find 


r=p+0 
r=p-0 


is  a  constant. 


(9)  W(m,k)  =  0  , 


A  =  [w(2m,2k)  ]"^j-(kad) 
3   2 


+  2  U,-(m^+k^d^)  I^(i-  r+  4  T„,) 


J  ^  "  ^2(j)  '"'  "^  "  '  ^({)^2    (f)   "   2(})' 
-  J  c^  -  2mc(mc+kbd) (K   +  2K      ) 

+  (mc+kbd)2[3  +  2^2^-  (m^kV)  K^(^%+  4^^^)]}  , 


] 


and  the  plasma  pressure  is  given  by 
(10)  p  =  i  (-a^+b^+c^)  +  J   e^  -k^a^d^ 


-1  +  (m^+k^d^)  ll 


+  3c^  +  4mc(mc+kbd)K .  +  (mc+kbd) ^ 


-l+(m^+k^d^)K^  1+  O(e^)  . 


*J) 


Finally,  a,  b,  c,  and  d   are  related  to  A,  B,  C,  and  D  by 


(11)  A 

(12)  B 

(13)  c 


aTTd^   +    e^aTTd^(^  -   k^d^    I    )    +   0(e-^)     , 


bTT(e^-d^)    +    e^TTd^i-   ^  b+ (mc+bkd)  kd   K 


+   0(e-^)     , 
Led    log   J-  +    e^Ld[i  c   -   ^  mkad   I.    +   i  m(mc+bkd)K 


l¥ 


<f) 


+0(e-')  , 


(14)  D  =  Lud^fl  +  I  e^  +  O(e^) 

We  have  calculated  the  equilibrium  up  to  second  order  in  e, 

as  the  stability  is  determined  in  that  order. 
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Ill .   Bifurcation  points 

Now  suppose  that  the  magnetic  fluxes  and  the  amount  of  plasma, 
(viz.,  A,  B,  C,  and  D)   are  prescribed  in  addition  to  the  size 
of  the  cylindrical  wall  (viz.,e  and  L).   Circular  configuration 
is  not  the  only  possible  equilibrium.   Helical  configuration  is 
another  possibility.   If  the  configuration  is  helical,  the 
equilibrium  is  as  described  in  Section  II.    The  parameters 
a,  b,  c,  d,  and  e  are  determined  in  terms  of  A,  B,  C,  D,  e,  and 
L  by  Eqs.  (9),  (11),  (12),  (13)  and  (14).   If  the  configura- 
tion is  circular,  the  plasma  column  has  a  circular  cross 
section   p  =  d',  the  magnetic  field  is  B  =  (0,0,a')   in  the 
plasma  region  and   B  =  (0 , c 'd '/r  '  ,b ' )  in  the  vacuum  region. 

The  parameters   a',  b",  c',  and  d'   are  determined  in  terms 

2  2    2 

of   A,  B,  C,  D,  e,  and  L   by   a'Ttd'   =A,   b'iT(e-d')=B. 

2 
Lc'd'log(e/d')  =   C,    and   Lird '      =   D.      Thus    the   two   sets   of 

parameters  are  related  by 


(15)  a'  =  a(l  -  e^k^d^  I   +  0(e-^)  |  , 

(16)  b'  =  b  +  e^(mc+kbd)k  d^  (e^-d^)  "-"-K   +  0(£^)  , 

(17)  c'  =  c  +  e^  -  ^  c  +  ^[-mkad  7^+c+  m  (mc+kbd)  K  |  (log  |)  "^ 

-HD(e^) 

(18)  d-  =  d[l  +  1  e^  +  0U'^)\  . 

When  the  circular  equilibrium  is  subjected  to  a  disturbance 
in  the  mode  with  wavenumbers  m  and  k,  its  stability  is 
determined  by  the  algebraic  sign  of 

(19)  (ka'd')^  I(d';m,k)  -  c'^  -  (mc'+kb 'd' ) ^K (d ' ;m,k)  . 
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At  the  bifurcation  point  (e  =  0) ,  a,  b,  c,  and  d  are  equal 
to  a',  b',  c",  and  d'   respectively,  and  expression  (19) 
becomes  W(m,k)   which  is  equal  to  zero   according  to  Eq.  (9) . 
Thus,  bifurcation  is  possible  only  when  the  values  of  A,  B,  C, 
and  D   are  such  that  the  circular  equilibrium  is  neutrally 
stable  for  some  mode  of  disturbance.   Since  near  the  bifurcation 
point  the  bifurcated  equilibrium  differs  from  the  bifurcating 
equilibrium  only  slightly,  the  bifurcated  equilibrium  is 
likely  unstable  against  the  very  disturbances  which  render  the 
original  bifurcating  equilibrium  unstable.   Therefore,  the 
bifurcated  equilibrium  can  possibly  be  stable  only  when  the 
bifurcation  point  is  in  the  neighborhood  of  the  margin  of 
stability  for  the  bifurcating  equilibrium.   The  circular  equili- 
brium is  at  the  margin  of  stability  when  W(m,k)  as  a  function 
of  m  and  k  attains  a  minimum  equal  to  zero.   It  is  known  '   that 
allowing  k  to  be  any  value  (ignoring  the  requirement  that  k  must 
be  a  multiple  of  2tt/L)  ,  the  minimum  of  W(m,k)  occurs  at  m  =  1 

or  m  =  0.   Furthermore  the  minimiom  of  W(0,k)  occurs  at  k  =  0 , 

2   2    2  2      2   2 
and  its  minimum  value  is   2a  -c  +2b  d   /  (e  -d  ) .   The  helical 

structure  disappears  completely  when  m  =  0  and  k  =  0  at  the 

same  time.   Therefore,  a  stable  helical  equilibrium  requires 

(20)  m  =  1  , 

and  the  parameters  must  be  such  that 

2 

(21)  2a^  -  c^  +  2b^  —^ — ^  >   0  . 

e  -d 

Also  the  plasma  pressure  has  to  be  positive.   This  requires 

(22)  -  a^  +  b^  +  c^  >  0  . 
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Let  us  consider  a  parameter  space  for  the  three  ratios  A/C, 

2 
B/C ,  and  Lue  /D   (or  equivalently  in  a  parameter  space  for  the 

2   2     2   2 

three  ratios   a  /c  ,   b  /c  ,   and   e/d,   see  Fig.  1  and  Appendix 

A) .   The  simultaneous  conditions  of  W  =  0  and  9W/8k  =  0 
(upon  elimination  of  k,  with  m  =  1)  define  a  surface  in  this 
parameter  space.   A  circular  equilibrium  represented  by  a 
point  in  this  surface  is  marginally  stable.   Circular  equilibria 
represented  by  points  in  the  region  above  the  surface  of  marginal 
stability  are  unstable.   Any  point  in  this  region  of  instability 
for  circular  equilibria   can  be  a  point  of  bifurcation.  However, 
only  a  set  of  discrete  points  will  have  the  associated  value 
of  k  equal  to  some  multiple  of  27t/L.   As  e  varies  among  a  class 
of  helical  equilibria  with  the  same  value  of  m  and  the  same 
value  of  k,  a  curve  is  mapped  out  in  the  parameter  space. 
The  point  with  e  =  0  is  the  bifurcation  point.   Such  a  curve 
may  stay  in  the  region  of  instability  or  may  penetrate  the 
surface  of  marginal  stability  and  enter  into  the  region  of 
stability.   This  depends  on  the  stability  of  the  corresponding 
circular  equilbrium  which  is  determined  by  the  algebraic  sign  of 

(a'-a)  ^+  (b'-b)  ^+  (c'-c)  ^+  (d--d)  ^. 

This  is  the  Taylor  expansion  of  (19)  at  the  bifurcation  point. 

As  was  mentioned  previously,  a  stable  helical  equilibrium 

must  have  its  bifurcation  point  in  the  neighborhood  of  the 

surface  of  marginal  stability  in  the  parameter  space.   Moreover, 

any  instability  for  such  a  bifurcated  helical  equilibrium  must 

arise  from  disturbances  with  an  azimuthal  wavenumber  m'  equal 

to  m,  viz.,  1  and  an  axial  wavenumber  k'  differing  from  k  at  most 

by  an  amount  first  order  in  e. 
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IV.   stability  Analysis 

According  to  the  energy  principle  for  stability,  an 
equilibrium  is  stable  if  all  disturbances  will  increase  its 
potential  energy.   For  a  hydromagnetic  configuration  with 

a  sharp  plasma-vacuum  interface,  the  second  variation  of  the 

3  8 
potential  energy  '   is  given  by 


where 


Wp  = 


^S  = 


0  <r<p 


6W  =  Wp  +  Wg  +  W^  , 


rp{v.|)^  +  $•$  1  dx  , 


r  (    1 

2 

r=p 


I     ^     Ip 


W   = 
V 


V$  .  Vl)  dr 


p<r<c 

t   is  the  displacement  of  the  plasma  due  to  the  disturbance, 
r  is  the  ratio  of  specific  heats,  and 


Q  E  V  X  (txB) 


n  is  a  unit  normal  vector  at  the  interface.   V<J>  is  the 
additional  magnetic  field  in  the  vacuum  region  due  to  the 
disturbance. 
Let 

N  =  (p,  -  3^  ,  -p  y|)    , 

which  is  an  outward  vector  normal  to  the  interface  so  that 
dS  =  N  de  dz.   Hence, 
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L    2Tr 


(23) 


W^    = 


0    0 


f       i    (N.t)' 
J       2  r=p 


->■ 
N 

•>    -»■ 
L  N.N 


B'B 


r=p+0 
Jr=p-0 


de    dz 


When 


satisfies  the  Euler-Lagrange  equation: 


(24) 


rp  V(V«C)  -  B  X  V  X  V  X  (^XB)  =  0  , 


Wp    is   minimized   to 

L     277 


(25)  Wp    = 


J    (rp   V-5    -   B.Q)^^       (N.t)^^      d0    dz    . 


0    0 


The  magnetic  field  V$   must  satisfy  the  boundary  conditions 


(26) 


^ 


r=e 


=  0  , 


(27) 


N«V$ 


=  N  •  V  X  (txB) 

r=p 


r=p 


Moreover,  since  V$   is  a  harmonic  vector  in  the  vacuum  region, 
$  satisfies  the  Laplace  equation 


(28) 


V^$  =  0  . 


Hence  W   can  be  written 


(29) 


W   = 
V 


L  277 

r 


0  0 


(_$)  _   (N.V$)  _   de  dz 


r=p 


r=p 


In  fact  W   is  minimized  when  $  satisfies  Eq.  (28)  . 

As  was  discussed  in  Section  III  we  only  have  to  scrutinize 
the  disturbances  with  wavenumbers  m'  and  k'  such  that 


(30) 


m'  =  1  ,     k'  =  k(l+es) 
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By  virtue  of  the  bifurcation  condition,  the  minimized  5W 
will  be  second  order  in  e.   The  two  sets  of  wavenumbers ,  one 
for  the  disturbance  and  the  other  for  the  helical  equilibrium, 
interact  and  generate  mixed  wavenumbers  appearing  as 
V(t>  +  v" 4)  '  ,  where  . 

(j) '  =   m'e  +  k"  z  , 

->■ 
for  C  satisfying  Eq.  (24) .   The  calculations  appear  differently 

for  s  =  0  and  for  s  7^  0 .     In  the  case  k*  =  k  the  displacement 

which  minimizes  6W  is  compressible,  whereas  in  the  case  k'  f^  k 

(more  precisely  when  s  is  zeroth  order  in  e)  the  required 

displacement  is  incompressible.   We  shall  consider  the  first 

case  in  this  section  and  the  second  case  in  the  next  section. 

It  is  unnecessary  to  scrutinize  the  gap   for   s    smaller  than 

zeroth  order  in  e,  if  the  ratio  d/L  is  first  order  in  e. 

Because  both  k  and  k'  must  be  a  multiple  of  2tt/L,  hence  their 

difference  must  be  zero  or  at  least  equal  to  l-n/L. 

To  simplify  notations  we  shall  put  d  =  1.   In  other  words, 
k  stands  for  kd  and  e  stands  for  e/d  henceforth. 

First,  we  shall  find  ^   which  satisfies  Eq.  (24) .   Since 
in  the  plasma  region  the  magnetic  field  is  given  by 

->■-»■       ->■      2->-      ,  3 
B  =  Bq  +  e  B^  +  e  B   +  0(e) 

with 

Bq  =  (0,0, a  )  , 

Bj_  =  a^(-iR,0,Z)exp(ict.)  ,    R  =  ^  T^  ,   0  =  E  I^  ,  Z  =  k  I^, 


^2 


=  J  ^2! -i  ar    h<p  '  ^    hi,  '  ^^  h<t>)  ^^p  (^24.)   , 
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we  shall  expand  ^    in    e   formally  as 


I  =  Iq  +  e  t-^  +  £  I2  +  0(e  )  • 


->   -t 


By  separation  of  variables,  E,        i      and  E,        can  be  resolved 
into  sums  of  terms  whose  dependence  on  G  and  z  has  the  form 
exp  (iv9  +  iKz) .   Only  certain  terms  of  such  form  are 
necessary  to  be  included  in  the  process  of  minimizing  6W. 
It  is  convenient  to  define  a  linear  partial  differential 
operator  L         (see  Appendix  B)  as 

(31)  L  X  E  rp  V(V«x)  -  Bq  X  V  X  V  X  (xxBq) 

->- 
Then  E,      must  satisfy 

(32)  '-  ^0  "  °^^^  • 

As  is  shown  in  Appendix  B,  various  harmonic  vectors  are 
orthogonal  in  contributions  to  5W.   Therefore  in  the  case 
k  =  k' ,  which  is  exclusively  discussed  in  the  rest  of  this 
section,  it  is  only  necessary  to  consider  the  particular 
solution 

(33)  Iq  =  a(R^^l  ,  ie«^^  ,  iZ^^J  exp  (ic^)  , 
where 

^0,1  ^   ^     U    '      "^0,1  "  F  ^(J)  '   ^0,1  =  '^  %  ' 

a  being  a  coirplex-valued  constant,  allowing  a  phase  difference 

between   B   and  !„ .   Since   V*?   =  0  and  VxVx(5^xb  )  =  0  , 

->- 

^,  must  satisfy  the  following  equation. 
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(34)  l|^  =  BqXVxVx  (IqXB^)  +  0(e)  . 

The  right  side  is  independent  of  0  and  z;  hence  so  is  the 
inhomogeneous  solution,  the  term  with  the  coefficient  a  in  (35) , 
Observe  that  among  various  terms  in  t^    only  those  containing 
exp  (14))  will  interact  with  f  which  contains  exp  {i<t>)    only 
and  contribute  to  6W  in  second  order  of  e.   Thus   C-,  must 
include  terms  containing  exp  (i2(|))  or  exp  (iO)  .       Indeed 
it  is  only  necessary  to  consider  the  following  particular 
solution. 

(35)  ti   =   J    (Ri^o    '°'°^    ■"    ^^4,0    '°'°) 

+  y(R][    2    '    ^®1    2     '    ^^1,2^     ^^P    ^^^*^     ' 
where 

-1,0    =   —2    ^-0,1^    -    ^0,1-)     ' 

1,0         2     ' 

1,2         dr       ^2<t>     '      ^1,2  F      '24)     '       ^1,2    ~    ^'^    ^24)     ' 

3  being  a  real-valued  constant  and  y  being  a  complex-valued 
constant.   Other  harmonic  vectors,  if  included  in  ^   will 
make  a  positive  contribution  to  5W  in  second  order  of  e, 
hence   may   be  excluded  in  minimizing  6W.   Now  5-j  must  satisfy 
the  following  equation. 

(36)  '-    ^2    =   BqXVxVx(|qXB2+|-lXB^) 

+   B^xVxVx(|^xb^+|^xBq)    +   0(e)     . 


-16- 


A  particular  solution  of 


L  X  =  B,xVxVx  (IqXB^+I^xBq) 


xs 


X  -  °' 

X   —   T 


(X^  ,-iXg  ,-iX^)  exp  {-i<t>)    +  J  (X^  ,iXQ  ,iX^)  exp  (i(j))  , 


provided  that  X   ,  Xq  and  X   satisfy 
^  r     G       z        -^ 


a\(^  X^-kX^)  +  (rp+a2)  d_  (1  d_  ^  ^^  _  m  ^^_^^^^ 


=      -a^k^(0S^  +  — ^  Z  S^)  , 
^    rp+a^     ^ 


m 


a\(^  X^-kXg)  +  (rp+a2)  ^  (^  ^^  r  X^-  '^  X^-kX^) 


2  2 

=   -  a  k  R  S„  , 


(rp+a^)k(-  ^  r  X   -  -  X.-kX  )  = 
^  r  dr     r    r   9    z 


2  2 
a  k  R  S   , 
z  ' 


where 


S„  =  - 


1  d 


r  3x 


dF(^0,l®  +  ®0,1^)J   ' 


z   dr 


i^  rCR^  ,Z  +  zJi 
r  or     0,1    0 


,iw) 


Hence 


^r^'^e-^^'^^^^^    ' 


X„  =  R  S„  +  -  S 

0       6    r   r 


X   =  R  S   +  k  S   , 
z      z      r  ' 


where 
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s   = 

r 


r    ,  I  (kr)K  (kT)  -  K  (kr)I  (kr) 

1   m     m        mm 

n  ^   l'(kT)K  (kx)  -  K'(kT)I  (kx) 
U    m     m        mm 


T(x)  dx 


is  a  regular  solution  of 


<F^-^-4-^')s.  =  ^'-> 


■9^S, 


On  the  other  hand,  a  particular  solution  of 


->■   ->- 


L    y    =    BqXVxVxC^^xB^  +  C-lXB^) 


IS 


=  T  (^r  ^  ^r  '  -^^e-i^e  '  -i^z-^^z)  ^""P  ^-^*^ 


^[K^r'i^e'i^z)^^(^2,l'iQll'i^ll^^  i(^2,l'i®2,l'i^2,l) 


exp(i<t)) 


(R^  .,10"  .,iZ^   )  +  ^(r^    i0Y    izY   ) 


a/  a 
T^^2,3 


2,3'""2,3 
provided  that  y  satisfies 


2  ^  2  ,  3  '   2  ,  3  ' -^  2  ,  3> 


exp  (i3({))  , 


-  a  1^^  ?  =  V  X  (IqxS^  +  l^xS^)  . 


Hence 


Y   =  R^  „(-  0  +  kZ)  , 
r     1 ,0  ^r         ' 


"^e  "  aF  ^1,0®  ' 


Y   =  - 


-  ^  r  R?  ,Z 
r  dr     1,0 


Z   =  - 


"^  a   ,2m     01,2.  -r    ,  -m  ^a    ,  ,  „a   >  d    t 

.^0,1^-7  -^   2k   )    I^^   +   (-0Q^^   +   kZg^^)    g^    I^^ 


^2 


z„  =  - 


2  d   ,  a    2m         a    d    ,   , 


Z   =  - 


2  1  d    -„a 

iT^  F  aF  ^(^ 


0,1  2^  ^2^  "^  ^0,1  aF   ^2(j)^  ' 
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where 


4,1  =  4,0^F^   ^^'^     ' 


®2,1  "  dF  ^1,0®   ' 


3   _  1  d     R 
^2,1  ~  F  dF  "^  ^1,0^   ' 


R 


Y 


Rjf  o(-  9  +  kZ)  +  (5-  Q^  _  +  kZ^  „)R  , 


2,1    ''l,2r 


^r  "1,2    ^^"1,2' 


e 


Y 


jY 


(R'  „0  +  0'   R  )  , 


Y 


2,1  -  3F  ^"1,2^  "  "1,2 


,Y 


^  "^   rCR^  ^Z  +  zT  oR) 


2,1    r  or     1,2      1,2 


Any  harmonic  vectors,  if  included  in  \^   as  a  homogeneous  part, 
will  mar-.e  contribution  to  6W  in  third  order  of  e.   Therefore, 
it  is  only  necessary  to  consider  the  following  particular 
solution 


(37)   ^2  =  ^  +  y 


where 


R 


\    (R2,-l'iQ2,-l'i^2,-l)  ^""P  ^-^*^ 


a  ,„a 


,a 


|(^2,l'iQ2,l'i22,l)  ^  ^(^2\l'iQ2,l'i^2,l"^ 


+  J(R5,i.iel,i,iz^,i) 


exp  (icf)) 


5"^^2,3'^®2,3'^^2,3^  +  2"(R2,3'^®2,3'^^2,3^ 


exp(i3()))  , 


a 
2,-1 

a 


=  x^  +  Y^  +  z^,  9^^_^=  -^e-^e-^e'  ^2,-1  =  "V^z'^z  ' 


R^  ,   =  X   +  Y   , 
2,1      r     r  ' 


0 


a 


,a 


2,1 


^e  +  ^e  '  ^2,1  =  ^z  +  ^z 
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w 


ith  the  above  solution  for  t,   we  have 


(38)   V'?  = 


,a  1  d    _,a   ,  „,     2  a  +a 


Tp+a 


^  ^S^R  exp(i(t))+0(e^)  , 


and 


(39) 


Q  =  aak(iR2^^,-02^^,-Z^^^)  exp  (i((>) 

+    g(0,0,-a)  +  2Yak(iR^  2'~®1  2'"^1  2^  ^^P  (i2(t')| 


+  e^  ^^-J^  ak(iX^,-X  , -X^+  — ^—^   S^r)  exp  (icj))  +  0(e"^)  . 

Tp+a 

Next  we  shall  find  $  which  satisfies  Eqs .  (26) ,  (27)  and  (28) 
In  the  vacuum  region  the  magnetic  field  is  given  by  Eq.  (8) , 
thus  a  straightforward  computation  gives 


(40)   N«Vx(|"xB) 


^^   =  iaEg^^  exp  (it}))  +  ie  (|  E^  ^2+'Ye{  ^3)  ^^P^^^*^ 


+  ie^f  J  E^^_^exp(-i(l))  +  (|  E^^^+BE^^^+  I  E^^^)exp(i<})) 


^  (J  ^2,3^  jE^^3)exp(i3<^) 


+  0(e-')  , 


The  various  coefficients,  E's  in  (40)  as  well  as  F's  in  (41) , 
A's  in  (43) ,  B's  in  (44) ,  and  D's  in  (46) ,  are  given  in 
Appendix  C.   In  view  of  (40),  4>  satisfying  Eqs.  (26),  (27) 
and  (28)  is 
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+  ie^ 


(41)   $  =  iaFg  ^K^    expCicf))  +  ie(jF^^2+  'T^l    2)^2(f)^^P  ^^^*^ 

+  (t^2,3  ^  J^2,3)S<^  ^^P  (i2*)J  ^  0(^')  • 

Now  in  order  to  evaluate  the  integrals  (23) ,  (25) ,  and  (29) , 
the  following  are  computed. 


(42) 


N    „  B«B 

: r-  'V  — =— 


N'N 


r=p-o  =  ^0^2-  ^C^exp(i(!)) 

+  e^  [^  C^+   ^  C3exp(i2(|)) 


+  0{e^)     , 


where 


Cq  -  -  c   , 


C,  =  2a^k^  +  8c^  +  4mc(mc+kb)K ,  -  2 (mc+kb) ^  , 
1  0  ' 


2,  2 


2.,  2, 


C2  =  -  2a  k  [1  +  2(m^+k  )  I   -  m^  I^j  +2c^ (-10+m^+k  ) 

+  8mc(mc+kb)  (1-2K'  )  +  2  (mc+kb)  ^  [l+2  (m^+k^)  K  -m^/C^  )   / 

C-  =  -  a^k^  fl+2(m^+k^)  I^  -  m^  I?]  -  c^  (10  +  8m\  ^^+m^+k^) 
3  1  9         (pJ  2(p 

+  4mc(mc+kb)  |  2-2K^+K  „^-2  (m^+k^)  fC^K  „^  1 
f     4)   2  9  q)  2  cp  J 

+  (mc+kb)^[l+2(m^+k^)K  -m^K^] 
+  4A[a^k^+c^+  2mc(mc+kb)K2(jj  -  (mc+kb)  ^J  . 
(43)  (N»|)^^p=  aA2^^exp(i(j))+  e  |  aJ  q+6aJ  O^^l^l  2^^^1  2)^^^^^^*^ 

+  £^||A^^_^exp(-i(t))  +(Ja^  3_+3A^  ^+  J  A^  ^)exp(i(t)) 
+  (|a^^3  +  J  A^  3)  exp  (i3())))  +  O(e^)  , 
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(44)   (Tp  V'?-  B'Q)^_^  =  a  B^   exp(i(J)) 


+  e 


+     £ 


1^1,0^  K,o-   (7^1,2-^^1,2)  -P  (i2,) 


jB^^_^exp(-i<^)  +  (Jb«^^+  6B^^^+  jBY^^)exp(i(t,) 


+  (|b^^3  f  lB^2,3)^^P(i3*) 


+  OCe-")  , 


(45)   (S.Dj^p  =  J-*<,A,1   -  7°^«^0,1^0,l^^P(i2(^) 


1   *  ,a   ,,a   ,^Y   ^  •  1 


a   .a 


J  a    otAQ^^(A3_^Q+A[^2^  +  I  ^^^^O  ,  1^1 ,0 


+  2a3A2^^Aj^O  +  a%A2^^AY,^2 


exp(i(()) 


+  (jaaA«^^A«^2  +  «Y^S,W,2)  ^^P  ^i^*)) 


+  e 


T  a  a(AQ^^A2^^+  2  ^1^0^1,0^  2  ^1,2^1,2>> 


+  ^  aa(AQ^^A2^_;L  +  j   ^1^0^1,0) 

+  a%(A«^,A2^^1  +  A^^qAJ^o)  -  I  °^^(^0, 1^2,1+  ^^,2^,2) 


+  B6  A^   A^    +1  Y*Ya^   a^ 
+  ^^  ^1,0^1,0  "^  2  ^  ^^1,2^1,2 


+  const,  exp  (i2()))  +  const,  exp  (i4({))   +  0(e  )  , 
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,a 


C46)   (-«>)^_,  =  iaD"  ,exp(i(|))  +  ie 


l^l,0+(|°l,2-^YDl^2)^^P(i2c^) 


+  le 


T  D2,-l^^P(-i*^  +  (t  °2,1  +  ^^2,1  -^  J  °2,1^  ^^P^^*^ 


+  (|D^^3  +  1-D^^3)  exp(i3*) 


+  Oie^)     , 


and  (N«V$)      is  given  by  (40) . 

r=p 

Only  the  constant  terms  independent  of  9  and  z  in  the 
integrands  contribute  to  the  integrals  (23) ,  (25)  and  (29) 
Thus  we  obtain 


(47)       6w  =   11^  a*a   W(m,k) 


+    £ 


2    ttL 


J  a   aW        +   J  aaW        +   agW        +   ^  a   yW    ' 


+    23B   W^^    +    Y*Y    W^"^ 


+   0(£-')     , 


where 


^        =   ^0,1^2,1   +   ^1,0^1,0    +   ^1,2^1,2    +   ^2,1^2,1 

+    C     (2A°'      A*^         +    A°'      a"  +    A^'      A^       )+    C    a'^        (a"      +A^       ) 

^0^^0,1^2,1    ^   ^1,0^1,0    ^   ^1,2^1,2^       ^1^0,1^^1,0    ^1,2^ 

+    ^2^0,1^0,1   +    °0, 1^2,1    +    ^1,2^1,2    -^    °2, 1^0,1    ' 


W 


,aa^      a         a  a  a         +   a"      B*^         +    C    (2A°'      a"        +   A°'      A°'      ) 

^0,1^2,-1    ^    ^2,-1^0,1    ^   ^1,0^1,0    ^    ^0^^^0,1^2,-1^   ^1,0^1,0^ 


+  c  a"     a"'       +  c  a"     a"       -  d^     e"         -  d'^       e^ 

^   '-l^O, 1^1,0    ^    ^3^0,1^0,1         "0, 1^2,-1         ^2,-1^0,1    ' 


w 


aB 


^0,1^2,1    ^   ^1,0^1,0    ^    ^1,0^1,0         ^0,1^2,1 


.a      ^3 


3 
^1 
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,a      „3       ,  t:,cx      _3 


+     2C     fA^^       A^  +    A"^       a'"'        ^     +    r    A*^       A^  +    n""       F'"       +E"'       D' 

^^0^^0,1^2,1    ^   ^1,0^1,0^     ^    ^1^0,1^1,0         ^0,1*^2,1^^0,1^2,1    ' 


W«Y  =  A?!  .bX  .  +  A«  .B^ 


+  B?  .aT  .  +  B^  .aT 


0,1  2,1    "1,2^1,2    "1,2"1,2    "0,1  2,1 
+  ^CqCA^^^A^^^  +  ^1, 2^1,2^  +  ^1^0,1^1,2 
^  °0, 1^2,1  ^  °1,2^1,2  ^  ^1,2^1,2  "^  ^0,1^2,1  ' 


63  -  A^   B^    +  C  A^   A^ 
^    -  ^1,0^1,0  ^  '^0^1,0^1,0 


^  (rp  +  a^)  -  ^  c^ 


w 


■YY  =  aY   rY 


.Y   aY 


Y   fY 


^1,2^1,2  ^    Vl,2^1,2  ^  °i,2^1,2  =  W(2m,2k)  . 


Finally  we  minimize  (47)  by  varying  the  relative  amplitudes 
and  phases  among  the  three  harmonic  vectors  in  ^.   Writing 


*   1/2  *   1/2 

a  =  (a  a)  ^  (cos  \l)  +  i    sin  i|j)  ,   3  =  (a  a)  ^  x  ,   y  =  a(x  +iy  ) 

P  Y      I 


and   utilizing  Eq.     (9) ,    we   have 


Pr,         2    ttL      * 
oW  =e      -»-  a   a 


l„aa         1   W     W  1   W    'W    '     ^/l   ,,,aa        1   W      W     \  _, 

.    2wM(x^^    1   ^  COS    *)2^   WYY^.    1   ^)2,   „TYy2j^   o(.3)  . 

With  W^^  >  0  and  W^^  >  0,  the  minimum  of  (47)  is 


(48)   6W  .   =£^11^  a*a 
mm       8 


,,aa   1  W^^W^^    1  W"]^fl 
^   -  T   „36   "  T   „YY 


-   W 


, aa    1  W  W 


^"^^ 


W' 


+  0(e-') 


This  is  the  minimum  for  6W  when  the  disturbance  is  predominantly 
in  the  mode  with  m'  =1  and  k'  ==  k. 
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V.   stability  Analysis  (Continued) 

In  the  preceding  section  we  considered  the  case  when  the 
predominant  mode  in  the  disturbance  has  the  same  wave- 
numbers  as  those  of  the  helical  equilibrium.   In  this  section 
we  shall  consider  the  case  when  the  axial  wavenumber  of  the 
predominant  mode  is  different  slightly  from  that  of  the  helical 
equilibrium  [m*  =  1,  k'-k  =  esk,  s  -  0(e°)  7^  0  ].  We  shall  use  the 
following  notations . 

I^,(r)  =    I(r;m',k')  ,  '^^^r)    =K(r;m',k')  , 

I^,_^(r)E   (k'-k)^d^T(r;0,k'-k)  ,  K^,_^(r)  e  (k  • -k)  ^d^fC  (r  ;0  ,k  • -k)  , 

I^jj,+^(r)  =  I(r;m'+m,k'+k)  ,        '^  <},•+(;,  ^^^  -  ^^  (r  ;m' +m,k' +k)  . 

The  perturbation  scheme  will  not  be  straightforward  as  in 
the  preceding  section  because  the  small  parameter  e  will 
appear  in  such  expression  as   exp  i(k'-k)z.   The  fact  that 
the  derivative  of   exp(ieskz)   is  one  order  smaller  than 
the  function  itself  causes  some  intermixing  among  terms  of 
different  order  in  e.      Thus  it  is  necessary  to  consider 
a  particular  solution  of  Eq.  (32)  with  the  following  form 


where 


with 


R 


^0,1  ^  °'K,l'^®0,l'i^O,l)  expdc^') 


0,1    ar  >•  '   ^0,1   T^  ^cf)'  '   ^0,1  "  ^  ^<t>' 
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is  a  harmonic  vector,  and  t„  is  an  undetermined  vector 

constrained  by  ^  f _  =  0(e)  and  V(V«f-)  =  0(e).   Although  f. 
9z   0  0  0 

does  not  contribute  to  6W  in  lower  than  second  order  of  e , 

>■  ->• 

its  inclusion  in  £,„  is  essential  for  the  expansion  of  E, . 

->• 
In  fact,  f   is  to  be  so  determined  that  the  first  order 

equation  is  solvable. 

Since        Lh^    ^   =   0,         Lf^   =  0(e)  ,      and     Vxyxd^xBg)    =   0(e), 

Eq.     (24)    becomes 

eit^  =  _  LfQ  +  |BQxVxVx|(hQ  j^+f q)  x  (B*+B^)  |  +  0(e)  . 

1   ->-        ->    -»-* 
In  the  right  side  the  term  j   eB  xVxVx(h_  n ^B, )  ,  which  contains 

the  factor   exp(ieskz),   is  not  in  the  range  of  the  operator  eL 

Therefore,  solvability  requires  that 

-  lIq  +  |BQxVxVx(hQ  ^xB*)  =  O(e^)  . 

This  is  achieved  provided  that 

->■->■      -»■ 

where   h.  -  is  a  harmonic  vector  constrained  by  -^  n.f^  r,  =  0(e)  , 
and   g^  _   satisfies 


Hence 


where 


^0,0  =  2  <^^0,0'^®0,0'^^0,o)  expdr-i^)  , 
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In  minimizing  6W  it  is  only  necessary  (any  other  permissible 
harmonic  vector  will  produce  a  positive  contribution  to  6W  in 
second  order  of  e)  to  consider 

^0,0  =  ^(^0,0'i9o,0'i^O,o)  exp(ir-ic^)  , 
where 

Then,  5^  satisfies  the  follov/ing  equation. 

->■ 
Observe  that  among  various  terms  in  ^„   only  those  containing 

->- 
expCicf)')  or  exp  (i({) ' -i(|))  will  interact  with  5-   which  contains 

exp(i(f)')  and  exp(i(f) ' -i({))  and  contribute  to  6W  in  second  order 

of  e.   Thus,  £;,  must  include  terms  containing  c-^xp(i(J) '  +  i(j))  , 

exp(i(})'),  exp(i(J) '-icf))  or  exp(i(}) ' -i2(())  .   Indeed  it  is  only 

necessary  to  consider  the  following  particular  solution. 

(51)    ti  =  fi  +  gi^2  +  ^1,1  +  ^1,-1  +  ^1,2  +  ^1,-1  • 

-y  -> 

Here   h,  „  and  h,   ,  are  two  harmonic  vectors 
1 , ^       ± , —1 

^1,2  =  HRl,2'iQl,2'i^l,2)  exp(ir  +  i*)  , 

\,-r    ^(^1,-1'"-®!,-!'^^!,-!^  exp(i(J)'-i2(j))  , 
where 


-27- 


^1,2  -  dF  >'+<})  '   °1,2 E"  >•+({)   '  ^1,2    ^^  ^'^^S'+c},   ' 

T,5    -  d   ,         o^    -  "^'"^"^  7        7^         -    fk'-2kn 

The  vectors   9i  o  '  ^1  1  ^"   ^1  -1  satisfy 

hence 

5l,2  =  7    ("1,2  -^91,2  -"1,2)  '^xPdO'+i*)  ., 

?1,1  =  II  K,l  '"1,1  'i^l.l)  *  J^"?,!  '"?,!  .i^Ll)l-P<i*')' 

where 

^1,2  =  ^^(90,1^-^0,1"^  ^  1^  ^  (^0,19-^0,1^)    =  °(^)  ' 

1,2  "  Fm?  F  3F  ''(^0,1^^0,1^^  ^  Fm?  -^r^  ^^0,1^  ^0,1^^  "^^^ ' 


Z 


^1,1  =  ^^  (^0,0^-^0, 0^>  ^  ^(^0,0^-^0,0^^  ' 
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Qi,-i  =  ^(9^,0^-^2,0^)  +  ]^-^^  (^2,0^+^0,0^)  ' 

„a     _    k    1  d    ,  a     a         ,      k   i"'-2m  ,  a     „a    . 
^1,-1  -  ^"-^7^   F  SF  ''(^0,0^^^0,0^)  ^  F-^^IT  -^ —  (^0,o"  "0,0^)  ' 

4,-1  =  1^^^^^  (4, 09^^0,0^)  -^  ^(4,0^-^^0,0^)  ' 

4,-1  =  ^(4, 0^-4, 0^)  -^  k^^  (4,0^-4,0^ )  ' 

^1,-1  -  ^^^rr^  F  dF  ^(^0,0^+^0,0^)  +  F^m:  -^  (^o,o°"°o,o^'   • 

And  the  undetermined  vector  f ,  ,  constrained  by  tt—  fi  -   0(e)  and 

-1-  o  Z    1 

V(V«f,)  =  0(e) ,  is  to  be  determined  by  a  solvability  condition 
for  the  second  order  equation. 

Since   VxVx(f^xBQ+  j   eh^  ^xb*)  =  0,  VxVx(g^  2''^0^  J  ^0  i''B3^)=0, 

VxVx(g^  l^'^O^  I  ^O^'^l)  ^  0'  ^"^^  VxVx(g^  -I'^^O'^  ^^  ^O^^^l)  "  °' 
Eq.  (22)  becomes 
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eht^   =   -   zLf^   +  I  BqXVxVx|(|^+^^^2+  91,1+  ^1,-1+^1,2+^1,-1^ 

Its  solvability  requires  that 

-  L  f^  +  |SoXVxVx|5^^^xg*  +  (^^^_^  +  h^,-l)^^J  =  0(e2)  . 
This  is  achieved  provided  that  ^]^  =  ^i  q  satisfies 

Hence, 

^1,0    ^  I  8<^^l,0'^®l,0'^'l,o)-'    4^^1,0'^®!, 0'^'l,o) 

+   f(^J,0'^^l,0'i^l,o)l      exp(ir-i*)     , 
where 

^1,0  =  -M(^i,i+^i,-i)^  +(^i,i-^i,-i^H  ' 
^?,o  =  ^MK,i+^?,-i)^-(^i, 1-^1,-1)^1 

+  ^~^Ie  |(^?,i+^?,-i)Q-^(Q?,i-'i,-i)^I    ' 

„a  -lid      ^l/oOi       ■r.o'         N7     .     ^701        7a         \r| 

^1,0  =  =      FaF  ^r ^1,1+^1, -1^^  -  (2^,1-^1,-1^^]   ' 

^Lo  =  -M(Q?,i+9i,-i)'-^'li^'l-i^'l 

+    ^"'   ir     |(4,X.,-l5®-(4, 1-4,-1^^1     ' 

70  -lid      ^l,pB      ,p6        X7    ,     f^B        7B        nrI 

^1,0  =  ^       F  aF  ^  I  (^1,1+^1,-1^2  +    (Zi^i-Zi,_i)R|    . 
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^1,0  =  -'^(^1, -1^-4,-1^)  ^  ^"'  ^F  (4, -1^-4,-1^)  ' 

Now  ty    satisfies  the  following  equation  . 
(52)   L   t^    =   2-BQxVxVx|(5^^2+5l,0^5i^_l+hi^2+hi^_l)xS* 

-^  (^1,2-^^1,1+^1,0+^1,2)^^1 

+  (^0,0+^0,1+^0,0)^(^2+^2)]  +  °(^)  • 

The  right  side  is  a  sum  of  five  terms  containing  exp  Cicj) '+iv4))  , 
V  =  -3,-2,0,1,2.   Only  the  inhomogeneous  solution  corresponding 
to  the  term  with   exp(i(()')  ,  which  is  denoted  by  g^  -^   and  satisfies 

^  ^¥52,1  +  I  ^M(?1, 2+^1,2)^^1+  ^l,0^^ll  =  °'   ^(^-^2,1)  =  °' 
will  contribute  to  6W  in  second  order  of  e.   Thus, 

where 

^2,1'Q2,1'  ^2,1'  ^2,1'  ^ll   ^^^    4,1    ^^'^   ^(^)  ' 


^2 


,1   =   ^^    (^1,29+^1,2^)  +  '^(^l, 2^+^1,2^)    ' 


e 


2,1  -    '^(^l, 22-^1, 2^)  +  F-aF  (^1,2^+^1,2^) 


-31- 


6    _  k   1  d    ,6   7  76   R^  +  ^  "^'  (z^   0-0'^   Z) 
^2,1  -  k^  F  SF  ^^^1,0^^1,0^^  +  F^  —  ^^1,0"  "1,0^'  • 

The  displacement  t  just  described  is  such  that 

(54)  V't  =  O(e^)  , 

and 

!2  ) 

Cx||  (BjL+B^)  +-|-  (B2+B2)  +  0(e-')| 


=  ^  §¥  I^O,l^^Q^o+  ^(^,2^h^,-l)  -^  °(^'H   • 


Accordingly, 

(55)   Q  =  a  1^  lho,i+hQ^Q+  ^  (^1 ,2+^1 , -1^  ^  ^^^^^ 
=  aak'(iR^^^,-0^^^,-Z^^^)  exp  (i({.') 
+  3a(k'-k)  (iR^^Q,-9^^Q,-zg^Q)  exp(i(f)'-i(}.) 
+  eyaCk'+k)  (iR^  2'~*^1  2  '  "^1  2^  exp(  i(j) '  +  i(f)) 


+  e6a(k'-2k)  (iR^  _i,-9i  -I'-^i  _;l)  ^^P  ^^*  '  ~^^*^  "*"  °^^  ^ 


And 


-32- 


(56)  N.Vx(|xB)  l^^p^Q  =  iaE^^^  exp(i())') 

+  ie[(jE^  2"^'i'^l  2)  exp(i(f)'  +  i(f))  +  6eJ  _^  exp  (icj)  • -i2(j))  ] 

+  (5-E2^Q  +  f  E^^q)  exp(ic()'-i(f))  +...]  +  O(e^)  . 

The  various  coefficients,  E's  in  (56)  as  well  as  F's  in  (57) , 
A's  in  (58),  B's  in  (59),  and  D's  in  (60)  are  given  in 
Appendix  D.   The  corresponding  $  in  the  vacuum  region  is 

(57)  $  =  iaF^^^K^,  expdr)  +  i  (aF^  ^  q+BF^  ^  q)  K^  ,  _^  exp(  i(|)  • -i(()) 

+  ie|(^F5;^2M,2)S'+c^  exp  (ic^ -  +  1(1))+  5fJ  ^_^K^  .  _2^exp(  i(}.  • -i2(j)) 

+  i^'l(|F^,l+  |f^,i  +  ^F^,i)K^.  exp(ir)+  ...1  +  0(e3)  . 
Now  we  compute  the  following, 

(58)  (N'l)^^p  =  aAp  ^exp(i(})')+  e\{^A^    2'^^^\    2)  ^^P  ^i* '+i*) 


+  (j  A^  o^^'^l  o)^^P^^*'~^'^)+  '^■^1  _iexp(  icj) '-12(1)) 

+  e^lQ-   A^^^+  I  A^^-^  +  J  A^^^)exp(i(f)')  +  ...|  +  O(e^)  , 

(59)  (rpV-l-B'Q)^^   =  aBp  ^  exp(i(f)')  +  £  [  (j  B^  2'^'^^1  2)  ^^P^  ^* '"^^^^ 

+  (|  B^^  0"^^^1  o)^^P(^*'-i*)  +  ^^1  _iexp(i(j)'-i2()))| 

+  s2|(|b«^^.|b3^^.XbY^^)  exp(ir)  +  ...)  tO(c-^)  , 


2  f/a 

^2,1  -^  2-^2,1  -^  t 
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(60)  {-^)^=      =    iaDg  ^  exp(i(f)')  +  i  (aD^  Q+gD^  q)  exp  (i(}) ' -i(})) 
+  ie|(jD^  2"^^^l  2)^^P^i*'^^*^^  ^°1  _iexp(i({)'-i2(t))  ] 


-    -'Uf  ■>2,l+l°2,l-  J°l,l)-''P<i*''^  •••]  +  °'^  '  ' 


and   (N»V$)     is  given  by  (56)  . 
r=p     "^      ^ 

Finally  we  obtain 


(61)   6W  =  11^ 


2,1 


aa ,  1 


,a3_,  1 


a  aW(m'  ,k')+c  (i  a  aW  +  ±.  a  6W  '^+  i-  a  yW 


^ay 


+  6*3W^^+  y*yW^^  +  6*6W'^^)  +  O(e^)   , 


where 


^    =  ^0,1^2,1  ^  ^1,0^1,0  ■'   ^1,2^1,2  +  ^2,1^0,1 

"^  ^0  ^  ^^0  ,  1^2  , 1'^^l  ,  0^1 ,  0^1 ,  2^1 , 2 ' 

+  c  a"'   (a"  +A^   )  +  C  A^  A^ 
^  ^1^0,1^^1,0  ^1,2^    ^2^0,1^0,1 

+  °0, 1^2,1  "-    °1,2^1,2  +  "^2,1^0,1  "^  ^°0, 0^2,0  ' 

aB  _  ,a    3    ,  ,a    B      a    6    a    B 
^    -  ^0,1^2,1  +  ^1,0^1,0  +  ^1,0^1,0+^0,1^2,1 

+  2C  (A^   A^   +a"   A^   )  +  C  A^   A^ 
+  "^^0  ^^0,1^2, 1^1,0^1,0^  ^  ^1^0,1^1,0 

-   -0,1=2,1  -  -0,1°2',1  -  '-tA.O   *   4,0»^,0  - 


w 


ay 


=  a"'   b"^   +A^   B^   ^-B^'   A^   +B°'   a"^   +2C  (A°'   A^   +a"   A^   ) 
^0  ,  1^2  ,  1^1 ,  2*^1 ,  2+^1 ,  2^1 ,2+^0, 1^2  ,  r  ^^0  ^^0  ,  1^2  , 1+^1 ,  2^1  , 2 ' 


+  C  A^'   A^   +  d"   e"^   +  D^   E^   +  e"   D^   +  E^'   d"^ 

+  ^1^0,1^1,2+  °0,1^2,r  °1,2^1,2+  ^1,2^1,2+  ^0,1*^2,1  ' 
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^        e  -1  ^ 
w"^^  =  W(m'+m,k'+k)  , 

w"^*^  =  W(m'-2m,k'-2k)  . 

Since  k  is  chosen  near  the  margin  of  stability,  it  differs 
only  slightly  from  k„,  at  which  8W/8k  vanishes.   Let 
J   (k'+k)  -  kg  =  eks'.   The  value  of  [W(l  ,k ' ) -W(l ,kQ) ] 
-  [W(l,k)-W(l,k-)  ]  gives  W(m',k')  =  e^s'sk^  ^ 


+    O(e^) 


k=kQ 


With  W^^  >  0  and  W^^  >  0,  the  minimum  of  (59)  is 
(")   ^^min  =  ^   ^aa^W   -  ^  -^  -  ^ 


W'^'"     ^    W^'^ 


+  4s'sk2  -i-^l  +  O(e^)  • 

dk    ' 


h°' 


This  is  the  minimum  for  5W  when  the  disturbance  is  predominantly 

in  the  mode  with  m'  =  1  and  k'  -   k(l+es) . 

We  remark  that  the  ratio  of  specific  heats  appears  in  (4  8) 

but  not  in  (62) .   Indeed  (62)  does  not  reduce  to  (48)  when  s  is 

set  formally  equal  to  zero.   This  nonuniform  behavior  is 

expected  since  the  perturbation  schemes  required  are  quite 

different  in  the  two  cases.   Presumably  if  an  intermediate 

perturbation  scheme  is  used  to  treat  the  gap  for  s  greater 

than  zero  but  smaller  than  first  order  in  e,  a  continuous 

variation  in  6W  .   v/ill  emerge, 
mm  ^ 
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VI .   stable  Equilibria  with  Long  Helical  Wavelength 

The  helical  equilibrium  in  our  discussion  will  be  stable 

if  both  the  expressions  (48)  and  (62)  for  6W  .   are  positive. 
^  mm     ^ 

In  this  section  we  shall  evaluate  them  in  the  limit  of  small 
kd ,  viz.,  when  the  helical  structure  has  a  long  wavelength. 

This  occurs  when  the  axial  fields  are  much  greater  than  the 

2      2     2      2 
azimuthal  field,  viz.,  a   >>  c  ,   b   >>  c  . 

When  the  argument  is  small,  the  modified  Bessel  functions  become 

I^(x)  ~    2"''(v!)"^x''  ,      K^(x)  *  2''"^(v-l)  Ix'"'  . 

Thus 

W(v,k)  =  1  a^<^  -    c^    +   k   (vc+bK)^  %-^    ' 

For   small   k      and  m  =  1  we  have  approximately  (evaluated 

at  r  =  d  =  1) 

1=1     I    =  ^ 

2  4 

K      -    -   ^   "^-^      k'        -        —   e  +1 

(j)     ~T~  '        '^2^  2  "TT  ' 

^      e  -1       ^       e  -1 

T=l       T      =7  T      -■'" 

^4)'     '     4)'-<J)     '    ^({)'+(})~J' 

2  4 

fc-   -    e  +1     K      _      2       ^      _le+l 

e  -1      T   T      e  -1       T   Y        e  -1 

With  the  asymptotic  values  for  the  various  constants 

A's,  B's,  C's,  D's  and  E's   as  given  in  the  ends  of  Appendices 

C  and  D  of  reference  10,  we  have 

-  \   c2+4c(c+bk)  ^l±|iti  _  2(c+kb)2  e^^^^f 
■y      _ e  -1  (e  -1)  (e  -1) 

2a^k^  -  c^  +  2(c+bk)^  V^ 

e  -1 

4                6   4     2 
aa_  o-,2,  2^„^2  e  +7    ^  /  ,  v,i  \  e  +e  +15e  -1 
^u,u,_  2a  k  +2c  —, -  4c(c+bk)  — ^ ■, 

e^-1  (e^-1) (e  -1) 
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+    2(c+bk) 


2    e^+2e^+30e^-2e^+l 
(e^-l)^(e^-l) 


e^-1  (e^-1) (e^-1) 


-    2(c+bk) 


2    4e^-2e^-lle^+l 
(e    -1) 


+    A 


4a V-    2c2-    8c(c+bk)    ^^1^+    4(c+bk)2    2e^2e2-l 


"^T 


,    2    ,.2 
(e   -1) 


o    2,  2,     2     .     ,    ^,,,     e^+1   ^,    _L,,x2    3e^+4e^-3 
3a   k  +c    -4c(c+bk)    — +(c+bk)       ^ 3— 


W 


.aB 


e   -1 


sk 


2,  2 


2 
2    e    +3 


6a''k''+2c"    "2  '  ""    -8c(c+bk) 
e   -1 


4_L-3    2    „    2 
e   +3e    -2e 


(e"-l) 

e' 
T 


+   6(c+bk) 


2 

4  ? 

2    e   +4e^-l 

2 1 

{e    -1) 


if   k'=k, 


if   kVk, 


„aY         ^r^    2,  2    ^    2    T  ,    ,    ^,  ,  ,     e^+e^+1   ^^  ,    ^,,.2    9e^+9e   +lle   -5 
W    '    =    10a   k   -2c    -16c(c+bk)    — 3 +2(c+bk) 


e   -1 


— 7 5 ' 


W 


33  _, 


1    ,p    _^    2>  12 

2-    (Tp+a    )     -   5-  c 


s^(2a^k^    +    -^  b^k^)     , 
e   -1 


if      k-    =   k    , 
if      k'    7^  k    , 


W^Y    =    2a V    -   c^   +    2(c+bk)2   £!±i 

e^-1 


Therefore 


(63)       6W 


mm 


2    ttL      * 
e      —3-  a   a 


aa      1   W^^w"^ 


W — - 


W 


■YY 


-   lw«"| 


if      k'    =   k    , 


(64^       «^min 


2    ttL      * 
e      —75-  a    a 


W      -   -^  


2,  2..  2,  2    e   +1n 


+    8  s  '  s  ( a   k  +b   k      -^5 — -) 

e    -1   J 


if   k^k, 
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Clearly  if  (63)  is  positive,  then  (64)  can  be  made 
positive  by  proper  restriction  on  the  admissible  values  Tor 
k'  ,  which  must  be  a  multiple  of  2tt/L.   This  is  so  when  2tt/L 
is  greater  than  the  difference  between  the  two  roots  of  k 
for  Eq.  (9)  . 

Now  in  the  limit  of  small  k,  Eq.  (9)  becomes 

/  2^,  2  e^+l\,  2  ^  „,  ^  e^+1  ^   ^      2    2   _  „ 

(a  +b  —^ — )  k   +  2bc  — ^ k  +  c   — =  0  . 

e  -1  e  -1        e  -1 

Its  two  roots  are 

(65)   k  =  ^  e^l  T    [(e^-l)(e^-1^2gQ)]'^'    ^ 
^  (e^-l)6Q-2e^ 

where   6q  =  1  -  a^/b^   is  the  plasma  beta  [1-a^/ (b^+c^) +0 (e^) ] 

2      2 
when   b   >>  c  .   The  expression  (63)  contains  three  dimension- 

less  parameters.   They  may  be  chosen  as  ^^,    e,  and  c/b.   Since 

c/b  is  regarded  as  small,  there  are  effectively  only  two  free 

parameters.   For  each  set  of  values  for  3^  and  e,  there  are 

two  values  for  k,  corresponding  to  two  different  helical 

equilibria  (with  different  values  of  L) .   In  general,  numerical 

evaluation  is  necessary  to  determine  the  algebraic  sign  of  (63) 

However,  some  simplification  is  possible  when  B^,  is  equal  to 

zero  or  when  e  tends  to  unity.   When  3^,  =  1/  the  two  roots 

2  2 

for  k  are  -c/be   and  -c/b.   With  k  =  -c/be  ,  (63)  becomes 

2*      2222      42 
e  a  aTTL(e  -2)  /e  (e  -1)  (e  -3e  +4)  ,  which  is  positive  for  all 

e  >  1.    On  the  other  hand,  with  k  =  -c/b,  (63)  becomes 

e  a  aiTL(-3e  +7- |2e  -4  |) /2(e  -1),  which  is  positive  for 

4 
1  <  e   <  11/5.   When  e  ■>  1,  (63)  becomes 
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e^a*a7TL3Q(l-6Q)  "■'"(8-266q-|4-136q  |)/16(e-l)  ^,  which  is  positive 

for  0  <  3o  <  4/13.   Simplification  also  results  when  the  two 

2 
roots  of  k  coalesce  at  e  -1+23^.  =  0.   Then,  the  value  of  6W  .  , 

0  '  mm' 

2  *     J  (e^-1)  (e^+4e'^-e^+8)-|e^-3e'*+2e^-2| 

e  a   aTTL 2 5 Z 2 ' 

(e  +1)  (e°-e  +7e^+l) 

is  positive  for  e  >  1.13.   But  this  special  situation,  namely, 
when  the  bifurcation  point  is  exactly  on  the  surface  of  marginal 
stability,  does  not  have  physical  significance  in  the  limit  of 
small  k,  because  the  corresponding  plasma  beta  is  not  positive. 

In  conclusion,  there  are  stable  helical  equilibria.   The 
shaded  regions  in  Fig .  2    show  the  parameter  ranges  where 
the  bifurcated  equilibria  with  long  helical  wavelength  are 
stable.   The  corresponding  circular  equilbria  with  the  same 
magnetic  fluxes  and  the  same  plasma  volume  are  mostly  unstable, 
as  indicated  by  the  undotted  regions  in  Fig.  2. 
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Appendix  A 

In  determining  the  marginal  stability  for  a  straight  sheet 

6  7  9 
pinch  extensive  numerical  calculation   '  '    for  the  simultaneous 

algebraic  equations  W(m,k)  =  0  and  dW/dk  =  0  were   required. 

The  laborious  computations  seem  to  have  arisen  from  solving 

the  transcendental  equations  for  e/d  and  kd  in  terms  of  a/c 

and  b/c .   The  difficulties  can  be  obviated  by  instead  solving 

2   2       2   2 
the  quadratic  equations  for  a  /c   and  b  /c   in  terms  of  e/d 

and  kd.   Elimination  of  b  yields  [I  and  1^   below  stand  for 

T(d;l,k)  and  K(d;l,k)  respectively] 


,a,4,6-T  dK   i^  dl,2    „,a,2,2 


2Ik3+  k(2K+k^)(I  §-K^) 


+  4K^  +  (2fC  +  k  1^)^ 


which  possesses  only  one  positive  root  for  (a/c)  ,  viz., 


,|)2  .  ,-4„  a<  .,  ai,-2  ,,,3,  ,„,,  K  «,  (I  «  -K   «, 


+  2K 
Correspondingly , 


lV*U2I  +k«)(I  «  -.«) 


1/2) 


j  • 


b 
c 


.-,z--<«,-M, 


I 


ZV.    k(2I.  k«)(I  ^-K«) 


i/^-k-V 


The  algebra  simplifies  in  the  following  three  asymptotic  regimes 

(i)   e  ^  1:   since  K  ~  -(1+k^)  ""'- (e-1)  ""'■  we  have 

a  k  I  -  c     0  , 
c  +  bk   =  0  . 

Hence  a^/c^  x   k"^!"-'"   and   b/c  z    -k"""-. 
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2      2-1 
(ii)   k  ^  0  ••  since  I  :::  1  and  K   z    -(e  +1)  (e  -1)    we  have 

2  2     2  2   2      2-1 

a^k"^  -  c   +  (c+bk)^(e^+l)  (e^-1)  ^   ^^    0    , 

b}\?-    +  bk(c+bk)  (e^+1)  (e^-1)  "■•-   =   0 

Hence   a^/c^  :;  2k"^  (e^+1)  "■'■  and  b/c  a  -2k"-'- (e^+1)  "■"■ . 

(iii)   k  -><»  [with  e-1  =  ©(k"""-)  and  k^(e-l)  <<  cosh^k(e-l)] 
since   I  ~  |k|~   and  K  :z    -k~  coth  k(e-l)   we  have 


2,  -1 


a  |kl-c  +  (c+bk)  k"-^  coth  k(e-l)  ^    0  , 


a^  |k|-2c(c+bk)k  "'■  coth  k(e-l) 


+  (c+bk)^[k  ■'■coth  k{e-l) -(e-l)csch^k(e-l)]  «  0 


Hence 


2 

a 

c 

and 


1 

nr 


1- 


^__^_  __^ 11/2|2 

^^^"^^"^^^^^      k^e-1)    V^V^   cosh2k(e-l)i    ) 


sinh  k(e-l)  cosh  k(e-l) 


i4i+  '^^^-'^ 


b 
c 


1 
k 


,,  sinh^2Me-l)  L  ,  J,  ^ 


k"(e-l) 


cosh  k(e-l) 


1/2  ( 


2   2  2   2 

Figure  1  displays  a  /c   versus  b  /c   for  various  values  of  e/d 


with  |kd|  as  the  varying  parameter  in  each  of  these  curves.  The 

dotted  curve  connects  the  points  on  these  solid  curves  where 

2   2    2    2 
2a  -c  +2b  /(e  -1)   vanishes.   The  straight  lines 

2  2   2 

a  -(1-3  ) (b  +c  )  =0  indicate  the  value  of  the  plasma  beta. 

It  is  seen  that  g-  is  in  the  range  between  1  and  4.99  [both 
W(l,k)  and  dW/dk  vanish  when  a^/c^  =  1,  b^/c^  =  0,  e/d  =4.99 
and  Ikdl  =  0.3  53]  . 
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Appendix  B 

-v   _ 

We  shall  solve  the  following  equation  for  x  =  (x  ,x„,x  ) 

L  X  =  0  . 

The  linear  partial  differential  operator  L      is  defined  by 
Eq.  (31) .   From  its  three  component  equations 

oX 

,„  ^  2,  1   8   n  ^  ^   2.  3^       19   8     .    rt 

(^P+^  )  F  9F  ^'^  ■"  ^  ^rr  ^e  -  F  ^  9¥  ^z^  =  °  ' 

t3X 

rp  l^vx  ^  0  , 

we  find  that 

3^       9   3      _    9 
7r"3Y37^z     "9r^' 


9x 


9^^     199      _19 

rr  ^e~F^^^z     FW^ 
3z 


a2 

V«x  =  J   g  , 

rp+a 


where   g(r,6)   is  an  arbitrary  function  independent  of  z. 

Upon  elimination  of  x   and  x^  ,  we  obtain 

n2  9        „2 
^   9¥  ^z  -  ^  5  . 

Hence 

9  ^  u  ^  9^H 

_  x^  =  g  +  h  +  -^  , 

o  Z 

where  h(x,9)  is  a  harmonic  function  independent  of  z,  and 
H(r,9,z)   is  a  harmonic  function  containing  z.   Accordingly, 
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2                 2 
T  ^r  "  TF  ^^    2^   ' 

3^       1  3   ,,  _^  9^H- 

r-2-  ^e  =  F  3¥  (h  +  r;?^   • 

3z  3z 

Integrations  with  respect  to  z  twice  yield 

X  =  f  +  z(F,G,g-'rh)  +  J   z^Vh  +  VH  , 

where  ?(r,e),  F(r,0),  and  G(r,e)   are  independent  of  x,  but 
constrained  by 

2 

Vf  +  g  +  h  =  — ^-^   g  , 
rp+a 

r  3r      r  39^ 

-*■  2  2    —1 

in  order  to  satisfy  V -x  =  a  (Tp+a  )   g.   Now,  periodicity  in  z 

requires  that  F(r,e)  =  0,  G(r,e)  =  0,  g(r,9)  +  h(r,9)  =  0, 

and  h(r,0)  =  const. 

Therefore,  the  admissible  solution  of  the  Euler-Lagrange 

equation 

L  I  =  0 

can  be  written  as 


I   =    y      y      a      J%-  I,     ,    i  ?■  I.     '    ikij     exp(ime  +  ikz) 

+      I      ^    (-i   -  F(r)  ,    |-  F(r)  ,    G(r))     exp (im0) +3 (|    ,f (r)  ,g (r) )  , 


where  F(r) ,  G(r) ,  f (r)  and  g(r)  are  arbitrary  functions  of  r 
The  corresponding  $  satisfying  Eqs .  (26),  (27),  and  (28)  is 
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•I"    =    y      y      ia      ,  (~  +kb)K.    exp(ime  +  ikz) 


[(|)"^^    (|)"^| 


,d,m       ,e,ml-l  f/r-m,     ,e,rnl  ,.    „, 

(-)     -    (;t)     I         Ks-)    +    (-)       expdme) 


Accordingly, 


5W   = 


_    ttL    ( 


^ 


U    J^    %,k%,k  W(m,k)d 


-2 


+      I     m^6*3F(d)F(d)  (-  ^  +  m 
m7«0  ^      d 


2  2m,  -,2m 

c  e      +d 

,2  2m    ,2m 

d  e      -d 


+    B3r(rp+a^)d^    -  ^  c^d^ 


) 


It  is  seen  that  in  this  case  various  harmonic  vectors  are 
orthogonal  in  contributions  to  6W. 
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Appendix  C 

The  various  coefficients  in  (43) ,  (44) ,  (46) ,  (40) ,  and  (41) 
are  given  below. 


>a 


0,1 


=      1    , 


A 


1,2 


2(in^+k^)I,     , 


A 


a 


o    2,  2 

2a   k 


1,0 


K     f 


vOt 


rp+a 

a 
2,1 

r' 
'2,-1  2,-1 


2,1 


2    2 
R*^    ,    -    2m^I,+    (m^+k^)  (l+I,)    +   ^^-^    [l+ (m^+k^)  I^]     , 


,a 


■4-' 


^  +  m^I,-    (m^+k^)  (i-  +   t  I.  +  Ar  ) 


rp+a 
'1^2    ^(j)  ■-(})■ 


^    2a^k^     r  2^,2>t2-1 


rp+a 


^1,0 


1 


ll  =  1  +  J  ("^'+^')^  ' 


^1,2 


1    , 


A^ 
^2,1 


(m2+k2)(I^+4l2^)     ; 


B 


a 


0,1 


a    k    I       , 


B 


a 


1,2 


=      a 


2,  2 


k^[2    -    (m^+k^)l|]     , 


B 


B 


a 


1,0 


2,1 


-   a   k    (m  +k    )  I  ,       , 


i^k^    -l-(m^+k^)I      +    2(m^+k^)^I^    + 


o    2,  2 
2a   k 

rp+a^ 


1   +    (m^+k^) if  I 
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'2,-1 


a^k^^  -    2m'-IJ+(m"+k")  I^(^  -21^^) 


.2,2,  ,_2.,  2,  ,     .1 


2    2 

4)       (J)         2(1)  r         ^ ^2      4)1  '  '({)J 


-    H'^    +    '("'^^')^^2j|     ' 


rp+a 


B 


B 


B 


B 


6 
1,0 

6 
2,1 

Y 
1,2 

Y 
2,1 


-      rp+a       , 


2,  2t 

-   a   k    I^    , 


2    2 


2c|>    ' 
2,_2ro         .  ,2.,  2, 


a   k^[2    -   4(m^k^)I    I      J     ; 


D 


,a 


0,1 


"  ^0,lS     ' 


D 


a 


.a      „ 


1,2 


-^1,2^    ^   <a[-^   ^    2(m^k^)K^K24)]     ' 


D 


a 


1,0 


=       -    E 


a 


0,1 


^a 


2,1 


-^2,lS    ^^l,2h    ^    ^(m^k^JK^K^J 

+   E2^Jl-2m2fc2+(m^k2)(2fC^+K2)-4(m^k2)2K2K2j     , 


D 


a 


2,-1 


D 


D 


'2,1 

Y 
1,2 

,Y 
2,1 


-    A[l    +    (m2+k2)K2]j      , 


"   ^1,2'^2<})     ' 

-^l,lS    ^   ^1,2  h   ^    2(m^k2)K^/(2j        ; 


<{>        2       (f) 
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e"    ,      =      mc   +   kb    , 
u ,  1 


,a  .  ^    .     . X  ,_2.,_2, 


E^    o      =      -    4mc   +    (mc+kb) (m  +k    ) (41 .-2K, )     , 
1,2  (()         ({) 

^2    1      "      mcr6-4(in^+k^)I    ]  +  (mc+kb)  [m^(-?I    +4K    )  +  (m^+k^)  (-1+1    ) 

2    2    ( 

-2(m2+k2)2i^K^  +  R-^^-j    +   £^   j-2mcl^+ (mc+kb)  [l+ 

+(m2+k2)I^(I^-K^)]j       , 
^2    -1   "      -mcr3+4(m^+k^)K2.1    +    (mc+kb)  [-m^  ( I    +2K    ) 

+(m2+k2) (3    +    1    1^+2^2^)-    4 (m^+k^) 2k^K2^-   R^ ^ J 

2    2    (  ) 

+    ^^   ^-      2mc7^- (mc+kb) [l   +    (m^+k^) I , ( I , -K , ) 1 
rp+a^   I  *  ^  (J)      <i)      4.    Jj 

+    Ar2mc+(mc+kb)  (m^+k^)  (I    +4K2.)1     , 
^2    1      "     -   mc   +    (mc+kb)  [1   +   j   (m^+k^)  (I    -K    )]     , 
E^    2      =      2 (mc+kb)     , 
^2   \      "     ~   ^"^^   "•"    ('^^^^^  (m^+k^)  (lA-'^(f)+4l2A)     ; 


F^*  =  E^ 

*0,1  ^0,1       ' 

^1,2  =  ^1,2   -    2e2,i(-'+'^')%    ' 

^2,1  =  ^2,1-    2E«^2("^'+^')So,-^^0,lKV("^'^^')(l^S) 
+    4(m2+k2)2K^K2^]     , 

pB  _  -pB 

^2,1  "         2,1    ' 


~47- 


^1,2      -      '"1,2     ' 

In   the    limit   of   small   k       [a/c   and  b/c   are   0(k~    )]    we  have 

3      2^,     ,    ^^,,    e'^+e^+l        oi^4.uv\^   e^  {e^  +  2e^  +  3) 
-    i  c   +4c(c+bk)    y, 2(c+bk)       — » 3 

^  e^-1  (e^-1)  (e  -1)  . 

X    «    — 5 ' 

2a^k^   -    c^    +    2(c+bk)^    V^ 

e  -1 


A°' 
^0,1 

= 

1  , 

A^' 
^1,2 

w 

2    , 

A°' 

X 

1    2,  2 
2a   k 

rr^-i-a2 

A^^l      =      0{k2)     , 

A^   _j^    !=      -    2    +    3X    +   O(k^)     , 

A^  =      i 

A^  =1 

^1,2  -^    ' 


A^  «      3    • 

^2,1  -^    ' 


T^oi  ^        2,  2 

^0,1  ~  ^   '^         ' 

T^oi  _        2,  2 

^1,2  ~  ^   k      , 

nCt  2,  2 

^1,0  -  -    a   k      , 


-48- 


B2    _-L    -      a^k^(2-3A)     , 


4,0  =  rP+^^    ' 

^2,1  ~  -    3   k       / 

I.Y  5      2,  4 

^2,1  =  -    6    ^   ^      ' 


2 

Cj^  :«      2a^k^   +   8c^  -    4c(c+bk)    ^y^  -    2(c+bk)^    , 

e  -1 

2  4         2 

C2  «.     -    4    a^k^   -    18c^   +    8c(c+bk)    ^^^  -    4(c-HDk)^    ^   ^^^  "^    , 

e  -1  (e  -1) 

„             ^           0.2,2           2    Ve'^-lS   ^    _    ,    ^u,  N     Se^-e'^-lle^-l 
C^  s:      -    2a   k     -    c      — 3 +   2c(c+bk)    — ^ ^ 

e^-1  (e  -1)  (e  -1) 

-    2(c+bk)2   ^It^e^ 

4 
+   4x[a^k^   +c^  -c(c+bk)    ^^  -    (c+bk)^]    ; 

e  -1  -' 

D«  .       (c+bk)    4^    , 

'  e  -1 

Dl^     ^      -    2c   ^  +    (C4bk)     3e^e^5e^-l    ^ 
e^-1  (e  -1)  (e^-1) 

°1,0      "     "    (c+bk)     , 

D-    ^      .      2c   ^'-f-l  -    4(c+bk)    ^'(V^%-^)     , 

D^  ~      -    c   ^=^  ^    ^^^^^j^j    ^'(e^-e2    2) 


.X[2c^-    (Cbk)    ^£^]     , 
e  -1  (e  -1)       J 
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D 


D 


D 


2,1 

Y 
1,2 

,Y 
2,1 


_  c  4^  +  (c+bk)  (^'^-^^^yi) 


e  -1 


(e'--l) 


(c+bk) 


e^l 
7^1    ' 


2                                        4         2 
T      e   +1    ,_,,,,.     2e   +3e  -1 
-    2c   -^ +    2  (c+bk)    2 ^— 


e   -1 


(e-'-l) 


^0,1 

,a 
^1,2 

.a 
^2,1 

,a 
'2,-1 

^2,1 


=      c   +  bk    , 


-    4c   +    2 (c+bk) 


•5    2    T 
3e   -1 


e   -1 


2c  -    4  (c+bk)    -5 —    , 

e  -1 

4c  2,    4      2     . 

e-5    ,     „-     ,i.i\       e    (e-e-4 
-    c   —3 —  +    2  (c+bk)    — -■' 

e  -1  (e   -1)  (e"-l) 


(e  -e   -4)     ,     ,  r-i       /     ,1,1  \     5e  -1   T 
-!^ 3 ^  +    A    2c- (c+bk)    — 3 , 

-1)  (e  -1)  ^  e  -1    J 


-    c   +    (c+bk) 


2e^-l 

2    ,     ' 
e  -1 


^1,2 


=      2  (c+bk)     , 


^2,1     ~ 


9    2_ 
-    2c   +    2  (c+bk)    =y-^ 

e  -1 
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A^ 

^0 

,1 

^1 

2 

a" 

,0 

A^ 

,1 

Appendix   D 

The   various   coefficients    in    (58) ,     (59) ,     (60) ,     (56) ,    and    (57) 
are   given  below. 

=      1    , 

=       [2    +   k'(k'+k)  jl^,     , 

=      1   +    k'2    +    k-k   +    k^    +    2    |t-  -    2    ^   I^    +    (-l+k'2_    2    ^)I^, 

+    k^(|T-+k^)I^    +     [2k(k'+k)+    k^(k'^  +  2k'k)    +    (2+k^)     |t-JI^,I^    , 
=      s   k    , 

=      1    , 

-       (|^+k2)I^.    +     [2(1+    ^)    +    (k'+k)2jl^.^^    , 

=     1     ; 


b"  -       a^k'^T 


4. 

0 

4. 

1 

A^ 
^1 

,2 

4 

,1 

4 

-I 

B^    ^      =      a^k'^    +   a^k'kfl  -    (l+k'k)I       I    ]     , 


5*^ 
1,2 


^1,0    =    ^^^'^  -  ^^^'H^  +  (i+^'^^^^r^J  ' 

B^  =.       =2i^|2f_-|     j_     M-i-V»2' 


^2,1      =      a\'^[-l   +    d+k-'^)!^,]    +   a^k'k|-2(l+k'k)  (I^,  +  I^)+4T^,I^I, 


^1,0     ~     ^^  ^'^<^'-<i>   ' 
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B^ 


l^^   =    sa'k[k'(k--k)    -    k    I^._^I^]     , 
^^   =      a2(k'+k)[k'    -    k(2+k'k+k2)T^I^,^^]     , 


^2 


B,       -,  —      a    (k  —2k)     ill    Ti     7 
1,-1  (j)'-2(f) 


"0, 

,1 

°"o, 

,0 

^l 

,2 

D« 

,1 

-^2,1%.    -^^1,2[-1    ^    (2-.k.2.k.k)K^.f(^.^J 

+   E^^  Jl+(2+k'2+k'k)K^.+     (-l+k'2)K2,-.(2+k'^k-k)2K2.K^.^^], 

^2,1=      -^S.j^(2-K^._^)    +(c+bk)[-I^._^+(l+k'k)K^] 

nY         -      _    E^       K 

°2,1    =      -    ^2,1%-    ^   ^l,2l-l   ^    (2+k'2^k'k)K^.K^.^^]     , 

n^  -  v^  k' 

°1,-1~       '     ^1,-1*^(1)  •-2(j)         ' 


^2,1    =       c    +bk.      , 

E^   2   =     -    4c   +    (c+bk')  (l+k'^)T      +(c+bk)  j^{  3+2k' ^+k  •  k)  I^  , 
-     (2+kk'+k^)K^]     , 
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+(c+bk')     l+k'2_(l-k'2)l^,+k^I^{k'I^.+kI^) 
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Figure  1 


Solid  curves  (parameterized  with  |kd|)  represent 

the  surface  of  marginal  stability  for  a  straight 

2   2    2   2 
sheet  pinch  in  a  parameter  space  for  a  /c  ,  b  /c 

and  e/d .   The  region  above  the  surface  is  unstable 
for  circulc^r  equilibria.   The  dotted  curve  repre- 
sents the  intersection  of  the  surface  of  riarginal 

2   2    2   2   2    -1 
stability  with  the  surface  2a  -c  +2b  (e  /d  -1)    =  0 

Thin  lines  indicate  the  value  of  the  plasma  beta 

2    2   2 
6=1-  a^/(b^+c^) . 
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Figure  2 

Helical  equilibria  with  lorg  helical  wavelength 
(kd  <<  1)  are  stable  in  the  shaded  regions. 
The  corresponding  circular  equilibria  are 
unstable  in  the  undotted  regions. 
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